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Abstract 
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properties. We also derive trace asymptotics and heat kernel estimates. 
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1 Introduction 



Stochastic methods based on functional integration applied to the study of properties of pseudo- 
differential operators and related semigroups offer a powerful alternative to the techniques of analysis 
[3118]. Typical problems addressed include spectral properties of the operator, heat kernel estimates, 
L^-boundedness, ultracontractivity properties, and the decay of the eigenfunctions. 

In the paper [lOj we study analytic properties of evolution semigroups generated by fractional 
Schrodinger operators 

Ho, = (-A)°/2 + < a < 2 

with (fractional) Kato-class potentials V . As a 7^ 2 these operators generate non-Gaussian a-stable 
processes running under the potential V . These are Levy processes with paths having jump discontinu- 
ities. The well-known case a = 2 corresponds to standard Schrodinger operators generating Brownian 
motion in the presence of V (which is an Ito diffusion under extra conditions on the potential). The 
properties of Schrodinger operators and fractional Schrodinger operators in many aspects markedly 
differ. One sharp contrast appears in the decay properties of their ground state (first eigenfunction) . 
Specifically, the ground state of a Schrodinger operator with pinning potential V{x) — )• 00 as — )• 00 
decays (super)exponentially while it decays only polynomially in the case of fractional Schrodinger 
operators. This difference is due to the heavy tails of stable processes, as opposed to the Gaussian 
tail of Brownian motion. Another remarkable difference is that, roughly, the Schrodinger semigroups 
g-t(-A+y) g^j^g intrinsically ultracontractive for super-quadratically increasing potentials, while this 
property holds for fractional Schrodinger semigroups e~*(^~'^^ ^ already for potentials increasing 
faster than logarithmically. 

Our aim in the present paper is to focus on one single pseudo-differential operator and derive 
fine details on its spectrum and eigenfunctions by using a combination of functional integration and 
hands-on analytic methods. We will consider the operator 

and obtain various formulae and estimates on its kernel, eigenfunctions, and spectrum. The interest 
in this particular choice is twofold. One is that in order to further develop the more general theory it 
is important to have cases of reference with as detailed information as possible. The results we obtain 
below are indeed more refined than the general methods using either pseudo-differential calculus or 
functional integration provide, and we view this paper as complementary to the more general results 
in [ZlllS]. A second motivation is that there is much controversy in the physics literature (see, for 
instance, \\.2\ |5]) about claimed solutions of fractional Schrodinger equations. Due to non-locality of 
these operators such equations are more delicate than usual Schrodinger equations and an appropriate 
rigorous mathematical treatment is necessary. The operator we consider describes the massless (semi-) 
relativistic quantum harmonic oscillator studied in physics. 

Our main results are as follows. First we derive a functional integral representation which allows to 
define H as a self-adjoint operator. From the results of [9j it follows that the first eigenfunction (ground 
state) ifi of H is bounded both from below and above by with suitable prefactors, moreover, for 
all other eigenfunctions |(/7„(x)| ^ const <^i(x) holds. Due to the special choice of the potential, by 
using special functions we improve this result to a detailed asymptotic expansion of eigenfunctions, in 
particular, tighten the order of magnitude on the bounds (Theorem 13.101 below) . Secondly, in it 
was proven that the eigenfunctions are uniformly bounded for the case of the Cauchy process run in 
an interval only. We prove uniform boundedness of all eigenfunctions for H onM (Theorem I3.13P , and 
show that the set of zeroes of each eigenfunction is finite (Corollary 13. 12p . Also, we discuss the shape 
of the ground state (Theorem I3.14p . On the eigenvalues our main results are the precise asymptotic 
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expansions in Corollary 13.61 resulting from Theorem 13.51 which we prove that the eigenvalues are 
simple. We also obtain a spectral gap estimate (Corollary 13. 9p . and derive the trace asymptotics in 
Theorem 13.81 We give a heat kernel estimate in Theorem 13. 131 

The plan of the paper is the following. In Section 2 we derive a Feynman-Kac-type formula for 
a class of operators covering our case and then particularize to our chosen operator. We use the 
functional representation to define the operator as a self-adjoint operator. In Section 3 we show that 
the Fourier transform of the eigenfunctions satisfy the Airy equation under appropriate boundary 
conditions. This allows us to identify the spectrum of the fractional harmonic oscillator operator and 
derive some asymptotic formulae. Furthermore, here we present the main results as discussed above. 
In a short Section 4 we provide an appendix of the used facts on Airy functions. 

2 Functional integral representation 

Recall that the linear operator with domain i/"(M^) = {/ G L'^{R'^) : G ^^(M'*)}, < a < 2, 

d ^ 1, defined by (— A)°/2/(A;) = \k\°'f{k), is the fractional Laplacian of order a. It is essentially 
self-adjoint on C^(M'^), and its spectrum is Spec((-A)"/2) = SpeCess((-A)"/2) = [o,oo). 

Let y : R'^ — ^ M be a bounded Borel measurable function. We call (— A)"/^ -|- F, < a < 2, a 
fractional Schrddinger operator with potential V, where V acts as a multiplication operator. Since V 
is a bounded function, the operator (— A)"/^ + F is self-adjoint on Dom((— A)"/^) defined as a sum 
of two self-adjoint operators. Therefore Spec((— A)"^/^ + V) C [0, oo). 

Let {Qx,J^x,Px) be a probability space and {Xt)t^o a real valued symmetric a-stable process 
on it, with < a < 2. {Xt)t^o is a non-Gaussian Levy process, in particular it has independent 
and stationary increments. We use the notations and E^, respectively, for the distribution and 
the expected value of the process starting in x G M at time t = 0; for simplicity we do not indicate 
the measure in subscript (while we do when have any other measure or process). The characteristic 
function of {Xt)t^o is given by 

EO[e^«^*] =e-*l^l", eeM, t>0. (2.1) 

As a Levy process, {Xt)t^o has a version with paths in D([0, oo); M"^), the space of cadlag functions 
(i.e., right continuous functions with existing left limits). 

Recall that a subordinator {St)t^o on a given probability space {^ls,J^S,Ps) is an almost surely 
non-decreasing [0, oo)-valued Levy process starting at 0. An example is the (a/2)-stable subordinator 
{St)t^o uniquely determined by its Laplace transform 

[e-^^*] = e-*^"^' , t ^ 0, A ^ 0. (2.2) 

Consider standard Brownian motion {Btjt^Q on a given probability space {VLw ^ Pw) ■, where 
Pw is Wiener measure. Clearly, 

Epwte'^^*] =e"*'^''' eelf^, i>0. (2.3) 

It is a standard fact that any symmetric a-stable process {Xt)t^Q can be obtained as a random time 
change of Brownian motion where this random time process is an (a/2)-stable subordinator {St)f^o. 
It is convenient to consider the processes {Bt)t^Q and {St)t^o on two diff^erent probability spaces 
{^}w,J^W,Pw) and (O5, J^5,Ps). Then the process {Xt)t^o can be obtained in terms of subordinate 
Brownian motion with respect to the (a/2)-stable subordinator: 

Xt : Up^. X Qp^ B {uj,t) i — > 55t(r)(w) := Xt{uj,T). 
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This can also be seen by the composition of the characteristic exponent ()2.3p with the Laplace exponent 
(|2.2p which gives (12. Furthermore, P can then be identified as the image measure of this process 
on L>([0,oo);M'^) such that 

P^{XteA) = {P^xPl){Bs,eA) 

holds for ah Borel sets A C M*^. 

In what follows we take the case d = 1, a = 1, i.e., the operators \J ^(^) = so that we 
will consider the fractional Schrodinger operator 



H:=[-^\ +x\ (2.4) 



dx^ ) 

The symmetric 1-stable process {Xt)t^o is also known as Cauchy process whose one-dimensional dis- 
tributions are given explicitly by 

p'^iXt £ dy) = --^rn^ — ^^y' ^ e ^• 

vr -I- (x — yY 

Our main concern in this paper is to study the spectral properties of H by using functional integration 
methods. 

First we show how H relates with the Cauchy process. We have the following Feynman-Kac-type 
formula, which we state in d dimensions and a class of V containing our special case (see O |10] for 
more general pseudo-differential operators). 

Theorem 2.1. Let F E L°°(M"') and {Xt )t^Q be a d- dimensional Cauchy process. We have 

(/,e-<v^+^)<7)= / dxE^\j{Xo)9{Xt)e-foy(^s)ds] ^2.5) 



Proof. We divide the proof into four steps. 
(Step 1) Suppose V = 0. Our first claim is 



(/,e 



f{Xo)g{Xt 



(2.6) 



We regard the process {Xt)t^o as the composition of Brownian motion {Bt)t^o and the 1/2-stable 
subordinator {St)t^Q as explained above. Let Ex denote the spectral projection of the self-adjoint 
operator —A ^ 0. Then by using (j2.2p and the usual Feynman-Kac formula for e*'^ we have 



''^^d{f,Exg) 



'XSt 



e-^*^d(/,i?A5) 



^Ps 



d{f,Exg) 
(/,e-^*(-^)^;) 



E' 



Ps 



dxE%JfiBo)g{Bs,)] 



dxE^[fiXo)g{Xt)], 



thus dMI) follows. 

(Step 2) Let = to < *i < ■•■ < in, /o, fn G L'^{W^) and assume that fj G L"^ 



We claim that 



{tj-tj_i)V^j. 



f{Xo)\{f,{Xt 



for j = 1, 2, n — 1. 

(2.7) 
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For simplifying the notation put Sj = tj — tj-i, for any j = 1, n and 

j = 1, ...,n - 1, gn = fn- 
Notice that gj = /je~*J+^^~'^(7j+i. By (|2.6p the left hand side of (j2.7p can be represented as 



9. = n{ U e-^^^f, 



f{Xo)gi{Xs,) -- 
Using (I2.6p again, we obtain 

E"" [gj{Xs^)] = / p{sj,y-x)gj{y)dy= / p{sj,y - x)fj{y)e-'^+^'^gj+i{y)dy 



Ey [p{sj,Xo - x)/,(Xo)5,+i(^s,+J] dy 

p{s„y- x)f,{y)Ey [g,+i{Xs^^,)] dy = E^ [/,-(X,^)E^^. J] 
for J = 1, n — 1. The above equalities yield 



f{Xo)MXs,)x 



and (I2.7P follows by the Markov property of {Xt)t^o. 

(Step 3) Let now O^V £ CbiW^). We show ([23]) for such V. Since is self-adjoint, the Trotter 

product formula holds: 



Combined with (Step 2) it yields 



hm /, e 



-(t/n)v/=Ag-(Vn)y 
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f,e 



-t{^/^+V) 



lim 



Since each cadlag path s i— )■ a;(s) = ^^(a;) is continuous in s G [0,t] except for at most finite points 
we have X]^=i(*/"')^(^ij/n.) ~^ lo y{Xs)ds as n ^ oo in the sense of Riemann integral. Thus ()2.5[ 
follows for V G 



(Step 4) We make use of the argument in [161 Theorem 6.2] to complete the proof. Suppose that 
V G L°°(M'^) and let K = 4'{x/n){V * where K = n'^()){nx) with G C^(M'^) such that 
0^0^!,/ 4){x)dx = 1 and (/)(0) = 1. Then Vn ^ V almost everywhere and Vn are bounded 
and continuous. Let M denote the set of all x such that Vn{x) does not converge to V{x). Then 
the measure of {t G [0, oo) : Xt{uj) G A/"} is zero P^-almost surely and jlyVn{Xs)ds — )• j'^V {Xs)ds as 
n — )• oo P^-a.s. Thus 



dxE^ 



f{Xo)g{Xt)e-^oyi^s)ds 



as n — 7- oo. On the other hand, e *(V A+y„) ^ t(V A+y) -^^ strong sense as n — t- oo, since 
+ K + y on the domain Dom(V^). □ 
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We can use Theorem l2.1l to define \/^^+V as a self-adjoint operator. We define the Feynman-Kac 
semigroup 

{Ttf ){x) = W [e-^o ^ f (z l2(M'^), X € M'^. 

Theorem 2.2. Let V G L^^ci^'^) such that V ^ 0. Then {Tt : t 0} is a strongly continuous 
symmetric semigroup. In particular, there exists a self-adjoint operator K hounded from below such 
that e"*^ = Tt. 

K can be identified as tlie self-adjoint operator \/— A + F for ^ y G ^[^^(M'^). 
Proof. Since ^ ^ G Llod^"^) t^ere exist &y \ c\)^ > such that 



sup 



(2.8) 



Therefore we have 



\Ttff ^ 



I 



dxE^ 



dxE- [\fiXt)\' 



with some Ct > 0. Thus Tt is a bounded operator from L^(M'^) to Similarly as in Step 2 of 

the proof of Theorem 12.11 it is seen that the semigroup property TtTs = Tt+s holds for i, s ^ 0. 

To obtain strong continuity of Tt in t it suffices to show weak continuity. Let f,g£ C^(M'^). Then 
we have 



{f,Ttg) 



dxE 



x,0 

PwxPs 



f{Bo)g{Bs,)e-f>^''^r)'ir 



Since S'4(r) — )• as t — for each r G fipg, dominated convergence gives {f,Ttg) — {f,g). 

Finally, we check the symmetry property T^ = Tt. Let Bg = Bs{uj,t) = — ^^^(^^(tj). 

For every t £ we have then Bg = Bg with respect to (= denotes that the random variables 
are identically distributed). Hence 



U,Ttg) 



dxf{x)E^^^p\e-^l>^^^s.)d^g{Bs,) 



0,0 

PwxPs 



dxf{x)e-^'' y(^+BSr)drg(^^ + ^^J 



E' 



0,0 

PwxPs 



dxfix - Bs,)e--^o yi^+Bsr-Bs,)drg^^^ 



In the second equality we changed the variable x to x—Bs^ Since Bs^ = —Bs^ and Bs,. — Bst = Bs^-Sr, 
we have 



{f,Ttg) 



dxE 



x,0 

PwxPs 



f{Bs,)e-foy(Bs,-Sr)drg^^^ 



Moreover, as St — Sr = St-r for ^ r ^ t, we obtain 



if,Ttg) 



dxE 



x,0 

PwxPs 



f{Bs,)e~^oy(Bs,_r)'irg^^^ 



dxE,- 



x,0 

PwxPs 



/(i?5je-/o^(^^^H5(x) = (ri/,5). 



The existence of a self-adjoint operator K bounded from below such that Tt = e *^ follows now by 
the Hille-Yoshida theorem. This completes the proof. □ 
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3 Eigenvalues and eigenfunctions of H 

3.1 Basic regularity properties 

Prom now on we consider H defined by (j2.4p and the related Feynman-Kac semigroup. {Tt : t ^ 0} is 
given by an integral kernel, i.e., there exists u{t,x,y) such that 



(rj)(x) = / u{t, X, y)fiy)dy, xeRJeL^ 



Lemma 3.1. For every t > 0, the operators Tt are compact. 
Proof. [9], Lemma 1. 



□ 



Lemma 3.2. Let u{t,x,y) be the integral kernel of the Feynman-Kac semigroup {Tt : t ^ 0}. The 

following properties hold: 

1. for every i > the function u{t, •, •) is continuous, strictly positive, and bounded on M x M 

2. the semigroup is intrinsically ultracontractive, i.e., there exists C{t) > such that u(x,y,t) ^ 



C{t)(fi{x)(fi{y), for all t > and x,y GM., where (pi is the first eigenfunction of 



Proof. i9j. 



□ 



By Lemma [3.11 above the spectrum of H is purely discrete and there exists an orthonormal basis in 
L^(M) consisting of eigenfunctions ipn such that Ttc/p^ = e~^"^ipn, where < Ai < A2 ^ A3 ^ . . . — )• oo 
are the eigenvalues. Using the relation between the semigroup {Tt : t ^ 0} and its generator H we get 
in strong sense 



Hipnix) = lim = hm 

40 t 40 t 



V^nix) = -\n<Pnix), X £ R, 



which means that the functions {(Pn)n&N are also eigenfunctions of the Schrodinger operator H and 



<Pn{x) + x'^ipnix) = \n(Pn{x), X £ R. 



(3.1) 



We conclude this section by discussing some basic regularity properties of the eigenfunctions of H. 
Lemma 3.3. For every n = 1, 2, . . . we have (pn G L"^(M). 
Proof. For every x E M define 



f{x)=E- 



and denote by r := inf{s > : Xg ^ — 1} the first exit time of the process {Xt)t^o from the half-line 
(— 1,00). Then, for every x G M such that x > 2 we have 



fix) 



^ E 



s$ E' 



-S^(Xs+xfds 



-{a-l)2{lAr) 



+ e 



g{z)dz, 
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where g{z) is a density function of the random variable r. The expUcit formula for g was derived by 
Darling [4] (compare also (Sj [TT] ) . There exists a constant ci > such that for every x > 2 we have 
g{x) < cix^ ([3], p. 286). Thus 

/oo / /"OO \ 



whenever x > 2. The same argument for x < —2 shows that there is a constant C3 such that for every 
|x| > 2 

/(x) ^ cse" 



Finally, we get 

\fniy)\dy = e^" / \Tnpn{y)\dy e^" / / u{l,x,y)\ipnix)\dxdy 



\x\^2 J\x\>2 

^ e^" ( / \ipn{x)\dx + C3 [ e"(l^l"^)Vn(2;)|c^a; ) < 00, 

./|x|>2 y 

where the last inequality is a consequence of the fact that ipn G (K) • □ 

In fact, a stronger property is true, which we will need below. 
Lemma 3.4. We have that x^Lpn £ L^{^), for all n = 1, 2, ... 
Proof. For V{x) = x^ the estimate 

C 

ipi{x) ^ 



F(x)|x|'^+" 

holds for the first eigenfunction (ground state) '^i(x), see [9j, Theorem 1. Since the Feynman-Kac 
semigroup is intrinsically ultracontractive by Lemma 13.21 above. 

Wn{x)\ ^ % \x\ > 1 (3.2) 

X* 

follows. Then the proof of the claim is straightforward. □ 
3.2 Eigenvalues 

We will determine the functions and corresponding eigenvalues A„ starting from the relation (j3.ip . 
Denote the Fourier transform of eigenfunctions by 

yxAx) := ^n{x) = I e''='^n{z)dz. 

v2tt Jk 



Note that yx^ is well-defined since (pn £ (1^) by Lemma 13.31 Moreover, Lemma 13.41 implies that 
yx„ E C^(M), n = 1,2, ... By performing Fourier transform in (13.11) 



- y'Lix) + \x\yx„{x) = Xnyxnix) (3.3) 
is obtained. We are looking for yx G C^{R) D L^{R) satisfying 

-y'i{x) + \x\yx{x) = Xyx{x), A > 0. (3.4) 
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Notice that if a function yx{x) is a solution of (j3.4p . then the function y\{—x) is also a solution. Hence 
it suffices to consider equation (|3.4p only for x > and construct even and odd solutions on the whole 
real line. 

For X > equation (13. 4p takes the form 

y'{{x) - (x - \)y\{x) = 0. 

On substituting z = x — \, f\{z) = yx{x) the equation reduces to the Airy differential equation (14. ip 
below 

r;{z) - zh{z) = 0. 

(For a discussion of Airy functions see the Appendix.) The general solution of the above equation is 
thus obtained as 

yx{x) = ciAi(x — A) + C2Bi(x — A). 

The facts that y £ L^{M.) and the function Bi(z) tends to infinity when z — )• oo (see (j4.3p ) imply 
C2 = 0. Without loss of generality we can assume that ci = 1. 

Finally, to obtain an even function on the whole real line we put 

yxix) = Ai(|x| - A), x£ M. (3.5) 

We furthermore require that the right derivative at is zero 

lim y'^{x) = 
a;— >-0+ 



which is equivalent to 



For odd functions we have 



and we require 



or equivalently 



Ai'(-A) = 0. (3.6) 
y^{x) = sgn(x)Ai(|x| - A), x G M (3.7) 



lim yx{x) 

!>0+ 



Ai(-A) = 0. (3.8) 



Conditions ()3.6p and (13. Sp together with equation (13. 4p imply that the functions yx defined on the real 
line by 1^ and 1^ belong to C2(R). 

By the Parseval equality and the fact that (xAi^(x) — (Ai(x))^)' = Ai(x), which is an easy conse- 
quence of the Airy equation (j4.ip , we get 



oo 

2 , 



ipl{x)dx = / {^nY{x)dx = 2j A\{x + \nYdx = 2[{k^-\n)Y + \nk\^{-\n)) . 

) J — oo 

We have thus proved the following result. 
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Theorem 3.5. The eigenvalues for the problem i3.1\) are given by 

hk-i = -a'k A; = 1,2..., 
^2k = —o^k, k = 1,2. . . , 

where and a'^ denote the zeroes of the functions Ai and Ai' in decreasing order. They are all simple, 
the eigenf unctions ip2k-i{x) are even and (p2k{x) are odd. Furthermore, the Fourier transforms of the 
L'^ -normalized eigenf unctions are given by 



^^^1^1-"") n = 1,3,5,... 



V2A;;Ai(-A„) 

, X G 

sgn(x)Ai(|x| — 



V2Ai'(-A, 



n = 2,4,6, 



Making use of the asymptotic expansions and estimates for the zeroes of the Airy fmiction and its 
derivative [HI [6] yields 



Corollary 3.6. We have 



^2fc-l ~ 9 ( ^vr(4A; - 3) ) , A; ^ oo, 
A2fc ~ / ( |vr(4/c - 1) ) , k^oo. 



where 



o/. f 7 2 35 4 181223 . 18683371 c 91145884361 , , 

^ V 48 288 207360 1244160 191102976 J ' ^ ^ 

2/,/ 5 2 5 4 77125 108056875 o 162375596875 mi ,^.r.^ 

^ ^ V 48 36 82944 6967296 334430208 ' ^ ' 



Moreover, we have 



o \ 2/3 

A2fc-i ^ (^(4A:-1) , A: = 1,2 



? -"5 • • • 5 



|(4.-1))"%A.. < (|(4.-l))"'(l4arctan(^^-^||. . = 1,2, 
Remark 3.7. A numerical calculation for the first few eigenvalues gives 



Ai 


r\j 


1.01879297164747 


A2 




2.33810741045976 


A3 




3.24819758217983 


A4 




4.08794944413097 


As 




4.82009921117874 


Ae 




5.52055982809555 



Using the above asymptotic formulae we can derive a result on the asymptotic behaviour of the 
trace of the semigroup at zero. 
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Theorem 3.8. We have 



lim t3/2ye-An* = ^. 

71=1 



Proof. We divide the series into two components 



OO OO 



k=l k=l k=l k=l 

By (13. 9p we get that for every t £ [0, 1] 



^ g-{3/2^(fc-3/4)2/3t) ^ ^^^^ ^ g-ct^g-(3/2^{fc-3/4)2/3t) 
k=l k=l 



for some constant c > 0. Moreover, the function x i— t- e (^/'^'"(^ 3/4)2/3t) g^p^^^^y positive and non- 
increasing on [1,00). Thus 



00 , / , roo 



{3/2.(x-3/4)2/3i)^^ ^ ^ e-(3/8.)2/3t ^ / ^_(3/2.(.-3/4)2/3t)^^ _ 



'1 V Ji 

A substitution yields 

fOO , 1 /-OO 



„-(3/2^(x-3/4)2/3t) , _ 1 / 



e 

(3/87r)2/3t 

and we obtain limi_j.o+ t^^'^F{t) = 7r^^r(3/2) = Iii the same way we get limj_^o+ t^^'^G{t) = 
and this completes the proof. □ 

Using the estimates for the eigenvalues given in Corollary 13. 61 we obtain an estimate for the spectral 
gap. 



Corollary 3.9. We have 



A2-Ai^ (I) '^'(3^/3-1). 



3.3 Eigenfunctions 



Next we derive the asymptotic behaviour of the eigenfunctions. We use the notation p^, Qn as in ()4.4p 
below. 

Theorem 3.10. For every k = 1, . . . and N = 2,3, . . . we have 



as \z\ —7- 00. 



For every k = 1,2, . . . and N = 2,3, . . . we have 

^2kiz) = V2l^— — + ... + (-1) ___J + o (^^2]vT3 J ' as\z\^oo. 
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Proof. For /c = 1, 2, . . ., we have X2k~i = ~o-'k-' ^"^^ 

I 9 1 r 



Ai(u + a'^) cos du. 



Integration by parts (2A^ + 2 times) together with (j4.5p . (j4.7p and ()4.4p give 
/ k\{u + a'f^) COS zudu = (-l)''Ai(^^^(u + 

iV 

+E^(-ir'Ai'-"(«+«;)°° + fl«w 

s=l 
^ 1 

s=l ^ 



where 



RNiz) 



COS 



z 



■2N+2 



Ai(2^+i)(^ + 4) 



oo 2^ 
+ 

z 



2N+2 



P2JV+iK)Ai(4 



1 



'OO 

Ai(2^+2) (^u + a'k) cos ziidu. 



Using the asymptotic relations (j4.5p and ()4.7p together with formula (j4.4p we get 

^ \^2 (|P2iV+i(4)Ai(4)| + jJ|"|Ai(2^+2)(n + 4)Mn) = 

Notice that pi(2;) = 0, which completes the proof for this case. 
For A; = 1, 2, . . . we have 



y/2 

^2k{z) = -TTT^ — T / Ai{u + ak)smzudu. 
Ai [ak) Jo 



Similar arguments give 



Ai>fc) 
^/2 



^2kiz) 



N N . 



s=0 
N 



(3.11) 



+ i2jv(z) 



s=0 
AT 



31+1 — («^-) + ^^(^) 



s=l 



where 



Rn{z) = 5S(-l)^+^Ai(2^+2)(n + a.) 



+ 



-1) 



N+l i-oo 



Z 



■2N+3 



Ai(2^+3) _^ Q,^) COS zudu 



= ^2N+3 ^g2Ar+3(Qfc)Ai^(afc) + j Ai^^^^^^ (u + Ofc) COS zudii^ . 
Thus we get \R]\j{z)\ ^ C2fc,Ar|z|~^^~^. Finally, notice that q2{x) = 0. This completes the proof. □ 
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Theorem 3.11. The eigenf unctions <pn OLf^ analytic functions on M. Their Maclaurin expansions are 
given by 



- \/34AiK)i; (2™)! 



2 1 ^a.2„,K)(-ir 2m 



\/2 ^ ■W2m+l(afc)(-l)"' 2m+l 

where /c = 1, 2, . . . and 



Ai (ofc) (2m + 1)! 



Wn{x) = I Ai{u + x)u^du. 

Proof. For A; = 1, 2, ... we have 



/ 2 1 f°° 
^2k-i{x) = W — TT^TTT / Ai(u + 4)cosxn(fu 



N2m 



^ -<Ma'k)Jo ''^0 (2m)! 

Moreover, by (j4.5p it is seen that there exists a constant > such that |Ai(M + a'^)| < Cfce~3"^''^ for 
all u > 0. Hence 



^ /"oo 2m 2m ^ |™|2m /-co 

-(2;;5r<-'>"'^'<"+°^"'" = iAi(»+<)i„-d„ 

°° |^|2m r 



2m roo 



'''^ (2m)! 



m=0 



By putting (im = (|) ^ T ( ^"^+^ ) and making use of Stirling's formula T{x) ~ \/27re 
as X — )• 00 we obtain 

= (2m + l)(2m + 2) ( -) , ^ 



- (2- + l)(2- + 2)^-J eV3^_^J ^__J 
Thus Pubini's theorem applies to get 



m— >oo 
> 00. 



^2.-i(x) = ^/2^ 1^ V AiTTT / Ai(n + a',).2-d J 



m=o V -«fc(2m)! VAi«) 

Similar arguments give 



□ 
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Corollary 3.12. Every eigenfunction (/?„, has a finite number of zeroes. 

Proof. Using the asymptotic expansions in Theorem 13.101 it is easily seen that for every n there exists 
An > such that sup|^|>yi^ > 0. This means that all zeroes of (pn are in [— Since the 

function is analytic, the set of its zeroes is finite. □ 

Theorem 3.13. The eigenf unctions ipn are uniformly bounded. 

Proof. We begin with the case of odd eigenfunctions ipn, where n = 2k — 1, k = 1,2, . . .; the proof for 
the even eigenfunctions is similar. It suffices to consider x ^ 0. We have 



/ 2 1 f°° 
^n{x) = A — - . , , , / Ai{u + a 1^) COS zudu. 

Using the asymptotic formulae ()3.9p . ()4.8|) for — a'^ and Ai(a'^) respectively we get that 

^ 0(A;-i/6) (3.12) 



V^Ai(a;) 
We have 

Ai{u + a'j^) cos xudu = / Ai{u) cos x{u + a'j^)du + / Ai{u) cos x{u — a'j^)du. 

The Airy function is non-negative on [a'^ , oo] and thus the absolute value of the second integral is 
uniformly bounded by Ai{u)du < oo. The first integral is a sum of two integrals Ii{x) and l2{x) 
where 

/■-< / sin i'^v?/'^ + -)\ 
Ii{x) = / Ai{u) ^ ^' cos x{u + ak)du 



r / N 1 f ■ f'^ 3/2 cosx{u + a'u) , 



Using the asymptotic expansion for the Airy function ([T] 10.4.60 p. 448) we get 

sin (1^3/2 + f ) / 1 



as ti — >• oo. Thus there is a constant ci > such that 



/oo 
u-'^/^du. 
-a' 



The integral 2,/ttI2{x) can be rewritten as the sum l3{x) + /4(x), where 



/3(x) = j (li^/^ + x) sin 0^^^^ + ^ + x-u + xa'fc^ 



(ui/2 + x)ni/4 



/4(x) = / sin -u ' + — — xu — xa' 
J~a[ V3 4 

The term J3 is uniformly bounded by Lemma |4 . 1 1 with 



/(x, u) = (n^/^ + x) sin Q""^^^ + ^ + xu + xa'^,^ 
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and 

g{x,u) = {u^/^ + x)~\-^/\ 

To deal with the term I4 we need to consider several cases. Let (x — 1)^ > — a'^ or (x + 1)^ < —a'l and 
rewrite /4 in the form 



An application of Lemma 14.11 with 

/(x, u) = (n^/^ — x) sin ^^^^''^ + — — xu — xa';, 

and 

g{x,u) = _ 2.)-i^-i/4 



implies that 1^ is uniformly bounded for x > y^— a^ + 1 and /c = 2, 3, . . .. For the case — a'^ < (x — 1)^ ^ 
— a'^ ^ (x + 1)^ we have 

U{x) = J_^^^ {u' -x)sm^-n/ +--xu-xa,>i _ ,) 

/■-< . /2 3/2 vr A d-u 

The first integral above is uniformly bounded by the same argument as in the previous case. The 
absolute value of the second integral is bounded by 

/■-< ^ ^ _ 3 3/4 _ 3/4 1/4 

The last inequality follows from the fact that (x — 1)^ ^ — and the asymptotic expansion for a^. In 
the case —a'l < (x — 1)^ ^ (x + 1)^ < —a'^ we split up 1^ as 

/4(x) = [f\ \[ [\{u^/^-x)^in(\u'l^ + \-xu-xa'^-J^ 



— a' 



3/2 , vr dn 



By Lemma [4. II the first two integrals are uniformly bounded following by the same argument as before. 
The last integral is bounded by C2k^^^ ■ Finally, when (x — 1)'^ < —a'^ < (x + 1)^ < — a'^ we have 



Ii{x) = I sin I -u ' + — — xu — xa^ 



' 4 --^) ^1/4 

/ / 1 /2 N • /2 o /o vr , \ du 

+ / yu' - x) sm -n / + - - xu - xa J ^ 



The absolute value of the first integral is estimated as 



4 n^/^^y-a; u^^- 
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The last term can be uniformly bounded by one more application of Lemma |4.1[ Hence we obtain 

/•oo 

/ Ai(n + a'l^) cos xudu = 0{k^^^) 
Jo 

uniformly in x ^ 0. Together with (|3.12p this implies that the functions (pn are uniformly bounded. □ 

Theorem 3.14. The ground state ipi is decreasing on (0, oo). Moreover, there exist xi > xq > such 
that ifi is concave on [—xo,xq] and is convex on (— oo, — xi] and [xi,oo). 



Proof. For < x < y, let {R[^^)t^o, {R\^')t^o be two squared-Bessel processes of dimension 1 and 
with index v = —1/2, such that i?Q^^ = x^ and i?Q^^ = y^, and let {r]t)t^o be a 1/2-stable subordinator 
independent from R^^^ and R^'^\ Using the comparison theorem (see [iSj, Chapter IX, Theorem 3.7) 
we get rI^^ ^ Rf^ for all t ^ with probability 1. Thus ^ , for all t ^ 0. However, the 
process R^^ is X| starting form x^, and is starting form y^. Hence 



Solids 



< X < y. 



For every x G M we have 



u{t,x,y)dy 

oo 

^e~^"^iPn{x)^n{y)dy 



n=l 



oo „ 

Ve-^"Vn(a;) / iPn{y)dy. 

n=l 



This implies 



ipi{x) 



fn{y)dy 



hm e^"* E^ 

>oo 



and therefore the ground state is non- increasing on (0, oo) as a limit of non-increasing functions. In 
fact, ipi is strictly decreasing on (0, oo). This easily follows from the monotonicity of ipi proven above, 
and the fact that ipi is analytic on M. 
Concavity follows from the expression 



/ 2 1 f°° 
'^i(O) = -J / Ai(u + a'^^du. 

V ^H^iJ Jo 



We have Ai(M + a'i)u'^ > on (0,oo) thus (^'/(O) < 0. By continuity of 99" at it follows that there 
exists xo > such that ^'{{x) < for all x S [— xo,xo]. Moreover, we have 



V9i(x) 



-a'l Ai(o'i) Jo 



Ai{u + a\)u^ cos xudu. 



Putting f{u) = Ai{u + a'i)u^ and integrating by parts we obtain 



1 f°° 
... , , / PAiu + a\ )u'^ cos xudu 
Ai(a'i) Jo 



/{5)(0) 



+ 



/' '{u) sin xudu, 
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where f^^\0) = 20Ai(a;) and f^'^\u) = P(n)Ai(n + a[) + Q{u)Ai'{u + a[), where P{u) and Q{u) are 
polynomials. Applying (14. 5p and (14. 7p we get 



lim xVi(x) = 20, 

Ixl— s>oo V ~Ol 



> 



and this implies that there exists xi > xq such that (p'l is positive on (— c«,— xi] and [xi,oo). This 
completes the proof. □ 

Theorem 3.15. The integral kernel u{t, x, y) is jointly continuous on (0, oo) x M x M. Moreover, there 
exists a constant c > 1 such that 



c(l+a;4)(l+y4) 
for every t > 1 and x, y € M. 



^ u{t,x,y) < c 



{l + x^){l + y*y 



(3.13) 



Proof. For every to > 0, using the uniform boundedness of the eigenfunctions ipn given in Theorem 
13.131 we have 

oo oo oo 

\u{t,x,y)\ ^ J];e-^"Vn(x)||(^n(2/)| ^Mj];e-^"*^M^e-^"*« < oo, x,yeR,t>to, 



n=l 



n=l 



n=l 



for some constant M > 0. The boundedness of the last series easily follows from the asymptotic 
expansions for An given in Corollary 13.61 Then the continuity of u{t,x,y) follows from the continuity 

of (fn- 

For A; = 1, 2, . . . we have 



V2 



V>2k-l{x) = ^ 1 + 



Ai«) 



kJ JO 



Ai^'*) (n + a'j^) cos xudu | , 



where K\^^\t) = t^Ai(t) + 2Ai'(t). From the asymptotic expansions ([Ij 10.4.60 and 10.4.62) one can 
easily get that 

|Ai(4)(t)Kcit2, 

with some constant ci. Thus, using (|4.5p and (|4.7p . we get 

/ Ai^^\u + a'f^) cos xudu ^ / \Ai^^\t)\dt + Ai^^\t)dt ci t^dt + C2 csi-a'f^f 

Jo Ja'. Jo Ja'. 



Similarly, we have 



where 



Ai>fc) 
V2 



-'P2k 



1 /"°° 

(x) = — 7 / Ai^^Uu + Ok) sin xudu, 
x'^ Jo 



Ai^'*) (u + ttfe) sin xudu 



< C3(-afc)' 



Using the asymptotic formulae (|4.8p and (|4.9p we obtain that \x^'^2k-i{x)\ ^ C4,{k — 3/4)^^/^ and 
\x^^2k{x)\ ^ C4(A:— 1/2)^^/^. This yields, combined with the asymptotic expansion for ipi, the estimates 



^2k-i{x) 

(/7l(x) 



^ c^{k - 3/4)11/6 and 



^2k{X} 

ipi{x) 



^C5(A:- 1/2)11/6. 
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for the ratio of the eigenfunctions, for every a; £ M. We have 

oo 



n=2 



Lpi{x) 



For every t > the series on the right hand side is uniformly bounded by 



n=2 



-(A„-A2)t 



ipi{x) 



'Pliy) 



^ C5 j^exp(-C6(A. -3/4)2/3) ^ 

oo ^ 

+ ^exp(-cr(fc- 1/4)2/3)/ A; -i 



k=l 



11 ■ 
1\~ 



Using the integral test for convergence it is easy to see that the above series converge. Thus the 
expression 



ipi{x) 



¥^n{y) 



¥^iiy) 



n=2 

tends to zero uniformly in x, y € M as t — )• oo. This proves the estimates (j3.13p for t > to and x E R, 
y G M for some to > 1- 

Now let t G [l^^o]- Then we have 

\u{t,x,y) - e-^'^ipi{x)^i{y)\ ^ cse-^'^ipi{x)ipi{y) 

with a constant cg > 0. This provides an upper bound for u(t,x,y), for all x,y gW. Using the above 
estimates for the ratio ipn{x)/ipi{x), dominated convergence and the asymptotic expressions for the 
eigenfunctions in Theorem 13. 101 we obtain that for every t G [l,io] 



lim 



u{t, x,y) _ , ^ exp(4t) 



|x|,|yKoo ipi{x)(pi{y) 



exp(a'^to) 



> 0. 



k=l 



k=l 



The function e is comparable with a constant on [l,to]) whence 



u{t,x,y) ^ ce ^^^ipi{x)ipi{y) 



for t G [l,to]) \x\ > xq and \y\ > yo for suitable xo,yo > 0. 
Notice that for every a G [—yo,yo] and s G [l,to] 



u{t,x,y) 



lim 



IE 



e"fc^ (f2k-iia) 

1 V-^ 



> 0. 



Positivity of the limit is a consequence of the well-known general estimate u{t,x,y) ^ ctfi{x)(pi{y) 
derived from intrinsic ultracontractivity, with a constant q > 0. Thus there exist es,£a > 0, Xg^a > 
and a constant Cs^a > 1 such that 

1 ^ u{t,x,y) ^ ^ 

for every {t, \x\,y) G {s - es,s + Eg) x (xs,a,oo) x {a- ea,a + £a). The family 

{(s-e„s + e,) X (a-ea,a + ea)}(,,a)e[i,to]x[-j;o,j/o] 
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is an open cover of the compact set [l,to] x [— yoj^o]- Therefore there exists a finite subcover 

{(■5fe-£sfc,Sfe + £sJ X (o-k - £ak,ak + Sak)}k=l,2,...,n 

of the set [1, to] x [-2/0, Vo]- Putting c = max{csj^,a^ : k = 1, . . . ,n}, xi = msx{xs^,ak ■ k = l,...,n} 
we get that 

1 ^ u{t,x,y) ^ ^ 
c ^ ipi{x)(pi{y) ^ 

for every (t, \x\,y) G [1, to] x (xi, oo) x [— yo, 2/o]- Due to the symmetry of u{t, x, y) we get the analogous 
result for |y|) G [l,to] x [— xojXo] x [?/i,oo) Since u and (^i are continuous and strictly positive 
we get 

u{t,x,y) ^ ce~^'^*ipi{x)ipi{y) 

for {t,x,y) G [l,to] x [— max{xo, xi}, max{xo, xi}] x [— maxjyo) yi}) niax{?/o, This completes the 
proof. □ 



4 Appendix: Airy functions 

For the convenience of the reader we summarize some basic properties of Airy functions used in this 
paper. 

The Airy functions Ai(x) and Bi(x) are defined as two independent solutions of the Airy equation 



/' — xy = 0, X E 



(4.1) 



The equation can be easily reduced to the Bessel equation (for x ^ 0) and to the modified Bessel 
equation (x > 0). This allows to express the Airy functions in terms of Bessel functions and 
modified Bessel functions K^, in the following way: 



Ai{x) 



1 he 



TT V 3 



K 



1/3 



„3/2 



, x^Q, 
x>0. 



(4.2) 



Bi(x) 




J-l/3 [ -^{-xf' 



2 , 



1/3 [ 3^ 



3/2 



+ ^-1/3 i 



3/2 



X ^ 0, 
X > 0. 



Using the relation Ai"(x) = xAi(x) we get that the n-th derivative of Ai is given by 

Ai('')(x) = p„(x)Ai(x) + g„(x)Ai'(x), 
where Pn and qn are nth order polynomials defined by the recursive relations 

Pn+lix) = p'nix) + XQnix), 

qn+i{x) = Pn(x) + g^(x) 
and po{x) = 1, qo{x) = 0. Below we give formulae for pn and qn for n = 1, . . . , 10. 



(4.3) 



(4.4) 
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Pi{x) 


= 


Qi{x) 


= 1 


P2{X) 


= X 


Q2{x) 


= 




= 1 


qsix) 


= X 




= x^ 




= 2 




= Ax 


qbix) 


2 

= X^ 


Pe{x) 


= + 4 


<l%{x) 


= 6x 


P7{x) 


= 


q-j{x) 


= x^ + 10 


psix) 


= x'^ + 28x 


Qsix) 


= 12x2 


P9{x) 


= 16x3 + 28 


Q9{x) 


= x^ + 52x 


Pio{x) 


= x^ + lOOx^ 


qw{x) 


= 20x3 + 80 



We recall some asymptotic results related to Airy functions. The asymptotic behaviour of Ai for 
large arguments is given by (p!], 10. 4. 59 and 10.4.60) 

Ai(x) ^ ^^^^^^^e-i^'^', x^oo, (4.5) 

, s 1 sin (1x^/2 +7r/4) , r,,. , , 

Ai(-x) = -^^^_^_^ + 0(^-7/4)^ (4.6) 

The corresponding formula for Ai' is given by 10.4.61) 

Ai'(x) ^ -^-^/^e-§^^' 
Asymptotic formulae for Ai(a'^) and Ai'(afc) are ([Ij 10.4.96 and 10.4.97) 



Ai'(x) - __L_xVV§x«/^^ (4.7) 



Ai(4) ~ (-l)'=-ivr-V2 (^^^ (A:-3/4)-V6^ ^ ^ (4.8) 

Ai'(afe) ~ (-1)^-^-1/2 /j (A:-1/2)V6, k ^ oc. (4.9) 

Finally we prove a result we have used in the previous section. 

Lemma 4.1. For any open set D C M. and a > Q let u ^ f(x,u) be a continuous function for 
every fixed x £ D such that f{x,u)du is uniformly bounded for (x,y) £ D x [a, oo). Assume that 
(x, u) I— )• g{x, u) is a uniformly bounded function on D x [a, oo) such that for every x £ D and y > a 
the function u — )• g{x,u) has continuous derivative on [a,y), and its derivative has at most one zero 
in [a, y) . Then 

ry 

F{x,y) = f{x,u)g{x,u)du 



is a uniformly bounded function in D x [a, oo). 

Proof. Denote by b the zero of the function g(x, •) in the interval [a, y). If there are no zeroes of g{x, •) 
in this interval h can be arbitrarily chosen in the interval. Fix x £ D. Applying the second mean 
value theorem for integration to the intervals [a, b] and [b, y] apart it follows that there exist constants 
Ai e [a, b] and A2 G [b, y] such that F(x, y) is equal to 

9ix,a) / f{x,u)du + g{x,b) / f{x,u)du + g{x,b) / f{x,u)du + g{x,y) / f{x,u)du. 

J a J \i J b J \2 

Using the assumption that the functions g{x,u) and f{x,u)du are uniformly bounded, the result 
follows. □ 

Acknowledgments: It is a pleasure to thank T. Byczkowski, T. Kulczycki and A. Strohmaier for 
discussions. 
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